Global Bifurcation Destroying The Experimental Torus T 2 
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We show experimentally the scenario of a two-frequency torus T 2 breakdown, in which a global bi- 
furcation occurs due to the collision of a torus with an unstable periodic orbit, creating a heteroclinic 
saddle connection, followed by an intermittent behavior. 



Chaotic behavior has been extensively studied in sev- 
eral areas of physical sciences [l[ , in economy Q , in ecol- 
ogy Q and in applied engineering One important 
question toward the understanding of chaotic motion is 
the transition from a regular behavior to a chaotic one. 
This transition has been studied in conservative [f| and 
dissipative systems [l|, Q . 

Among the ways chaotic behavior appears are the 
routes via the bifurcations of a quasi-periodic attractor 
with N incommensurate frequencies, also known as torus 
T . This route for the higher dimension torus, in which 
N > 3, has been well analyzed in Refs. For the low 

dimension torus T 2 chaos might appear by the Curry- 
Yorke scenario [H, [TO, EH , and by torus doubling bifur- 
cations [Hj]. A torus T 2 transition to chaos was also 
observed in a COi laser fl3j 

Recently, Baptista and Caldas [3, EI| have shown the 
appearance of chaotic behavior by an abrupt transition 
directly from the torus T 2 , in which there is a collision of 
the torus with saddle points [l6j originating an intermit- 
tent behavior. They argued that there is a mechanism 
of reinjection associated to a heteroclinic connection be- 
tween the central unstable focus and the saddles, which 
takes the trajectory back to the vicinity of the focus. 
As reported by Letellier et al.[17j. this phenomenon is 
not only restricted to continuous-time systems but it also 
happens in spatio-temporal systems. 

In this work, we report experimental details of this 
scenario of T 2 breakdown. We show the global dynamics 
by identifying the collision of the torus with the saddle 
points, and by characterizing the local dynamics show- 
ing the existence of the focus and saddle points, i.e. the 
stretching and the folding character of the saddle vicinity, 
and the spiraling character of the focus vicinity. Finally, 
we also experimentally characterize the observed inter- 
mittency. 

The forced Chua's circuit is composed of two capac- 
itors, Ci and C2, two resistors, R and r, one induc- 
tor, L, and the non-linear active (piecewise-linear) resis- 
tor, _Rwr.[l4j]. whose characteristic curve is mathemat- 
ically represented by: inl{V c i) = moKi + 0.5(mi — 
m n ){\V c i + B p \ — \V C \ — B p \}, where V c \ is the voltage 
in the capacitor C\ . The forcing applied to the circuit is 
of the form V p sm(2nf p t), where V p is the amplitude and 
f p is the frequency. More details about this circuit can 



be found in Ref. [IS]. 

The circuit is constructed such that the char- 
acteristic curve presents mo=-0.539(3)mA/V , m\—- 
0.9lO(2)mA/V, and B p =l. 200(4) V. The number be- 
tween brackets indicate the deviation in the last digit. 
The components of the circuit are C\ — 0.0047/iF, Ci — 
0.052/xF, R € [1.0, 1.7}KQ, L = 9.2mH, and r = 10S1. We 
set the resistor R to obtain the regime of a Double-Scroll- 
like attractor, and then we introduce the perturbation 
destroying this attractor, which gives place to the attrac- 
tors observed in this work. We use a Tektronix AFG320 
function generator connected in parallel to the resistor r, 
to introduce the forcing in the circuit. Data acquisition is 
performed using an AT-MIO 16E1 National Instruments 
board (12 bits) with a sampling of (5=l/250kHz. 

We acquire the variable V Cl (t = iS), with i representing 
the number of acquired data points. Next, we identify 
all its local maxima V Cl (n) , where the n index represents 
the number of maxima. This time series V c \ (n) is used 
to reconstruct the attractor V Cl (n) x V Cl (n+p), using 
time-delay embedding coordinates, with p representing 
the time-delay for the embedding space. The value of 
p is chosen conveniently according to the geometry of 
the attractor studied, here p = 2. By working with this 
embedding, we can better visualize the attractors and the 
manifolds since this space has one dimension less than the 
usual phase space. 

We use as control parameters the voltage V p and fre- 
quency f p of the forcing. As we change these parameters 
several windows of periodic motion appear. In particu- 
lar, we explore the system close to the periodic window 
at V p = 0.7V and f p = 3500Hz. For V p = 0.7V fixed, we 
change f p within the interval [3505,3537] Hz and build 
the bifurcation diagram as shown in Fig. [1] The periodic 
window remains up to f p = 3514 Hz, when a Hopf bifur- 
cation occurs. Further increasing the control parameter, 
at f P — 3530.6Hz (indicated by letter C) a sudden change 
in the system dynamics happens, the destruction of the 
torus by a torus-saddle collision, generating an intermit- 
tent behavior. Close to the Hopf bifurcation point, the 
torus has a circular shape. As the parameter is increased 
the torus shape becomes like a pentagon. This change 
happens when the torus grows and gets close to the five 
fixed points corresponding to a period-5 saddle orbit. 

The existence of the saddle points, near this distorted 
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FIG. 1: Bifurcation diagram of the forced Chua's circuit, 
for V p — 0.7 V. One sees a fixed point that bifurcates by a 
Hopf bifurcation into a torus T 2 . At a critical parameter f p — 
3530.6Hz, indicated by the letter C, a torus-saddle collision. 
The attractors indicated by the letters A and B are shown in 
Figs. [2] and [3] respectively 



FIG. 2: (Color online) The pentagon-shape structure is the 
torus T2, and with the small filled circles we show orbits that 
are placed in the vicinity of the saddle as we perturb the 
circuit. This takes the trajectories away from the torus and 
reveals the structure of the phase space near the torus, i.e., 
the existence of the saddles. 



torus, can be verified by introducing a small perturbation 
in the forced Chua's circuit, e.g. changing the magnetic 
flow in the inductor. This perturbation sets the trajec- 
tory apart from the attractor[ll|. 

With this procedure we are able to place the trajectory 
in the vicinity of the saddle points, as well as, in differ- 
ent parts of the phase space. We repeat this perturbation 
until we have a picture of the phase space in which we 
can clearly see the saddle points. Those are identified 
by selecting trajectory points that after p-iterations are 
mapped close by. In Fig. ® we show the torus that re- 
sembles a pentagon and the perturbed trajectories that 
are placed in the vicinity of the saddles. After reaching 
the saddle vicinities, the trajectory goes either toward 
the torus or away from it, by the stable or the unstable 
manifolds of the saddle points (filled circles), respectively. 
At a critical parameter f p — 3530.6 Hz, the torus collides 
with the external saddle points, it breaks and an intermit- 
tent behavior takes place. Analyzing the destroyed torus, 
depicted in Fig. A) , we can still detect the presence of 
the saddle points, which indicates that at the collision of 
the saddle with the torus T 2 a heteroclinic saddle con- 
nection is formed. That is done by analyzing the local 
geometry of the mapping around the saddle points, i.e. 
observing the return of a set of initial conditions around 
these saddles, after 5 iterations. In Fig. \*2B), inside 
the square shown in the Fig. [3]^A) we connect with ar- 
rows the points in the set of initial conditions and the 
corresponding returning points. We selected the shortest 
arrows in order not to overwhelm the picture, revealing 



the saddle. 

We also see that the dynamics is still related to the 
torus. In fact, the trajectories behave as if they were in 
a quasi-periodic oscillation, near the saddle points, for 
a while. Then, they go away from the saddle vicinity 
along the saddle unstable manifolds, spending some time 
far away from the region delimited by the saddles. Fi- 
nally, they are eventually reinjected in the unstable focus 
vicinity, spiraling toward the saddle points. 

The unstable focus, in the attractor shown in Fig . 
|3JA), was obtained by the Fixed Point Transform [2CJ], 
which consists in transforming the original data set, such 
that the transformed data is concentrated on the pe- 
riodic orbits. The histograms of the transformed data 
presents sharp peaks at the location of the unstable pe- 
riodic orbits, which allow us to estimate their positions. 
The trajectories in the branches of the broken torus, the 
extended structures in Fig. [3JA), are governed by the 
unstable manifold of the saddle points. To see this, we 
evolve in time a set of points in the branch, close to the 
saddle, and after k x 5 iterations, with k small, this set re- 
turns to the branch, but stretched. Following the consec- 
utive iterates of the branches, we notice that the points 
on the branches are eventually mapped in the vicinity of 
the focus. The scaping of the trajectory from the saddles 
(chaotic bursts) and the subsequent reinjections around 
the focus generate the intermittent behavior. 

In experimental situations, measuring the average lam- 
inar time, namely r, can be rather problematic, as noted 
in Ref. [2l| . since one might mislead the determination 
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FIG. 3: (Color online) Fig Eft A) The chaotic attractor 
born from a global bifurcation, the torus-saddle collision, for 
f P = 3531Hz. The square in the center of the broken torus 
represents the unstable focus point. Fig. EftB), we analyze 
the dynamics inside the small square in Fig Eft A). The ar- 
rows connect a given point with its 5th iterate. One can see 
the saddle structure represented by the arrows, the unstable 
(stable) manifold is revealed by looking at the direction of 
the vectors. The filled circle indicates the saddle point posi- 
tion and the oversized arrows indicate the unstable and stable 
manifolds of the saddle point. 



of the laminar period, by working with variables that are 
in a projection of the full space. It is better to com- 
pute a more robust quantity called 7. This quantity 
measures the fraction of the time of the chaotic burst 
7 = tburst/Uotai and, as the chaotic burst is very pro- 
nounced, the error in the statistics is smaller. Assuming 
r oc e~^, the variable 7 is related to r by 7 = (l + ^) 
with d being a constant (see Ref. [21| for more de- 
tails). Since 1 << d/e^, we can write 7 oc e M . From 
our experimental results we have that fi = 0.96 (Fig. 
H|). Next, we estimate the Lyapunov exponents of our 
data through the Eckmann-Ruelle technique ((J, by us- 
ing the two-dimensional embedding {V c i(n), V c i(n + 2)}. 
The spectra of the Lyapunov exponents, for the forced 
Chua's circuit has the form A = (Ai, 0, A2, — c). The 
null exponent corresponds to the direction of the flow, 
and the — c exponent corresponds to the strong stable 
foliation. The other two exponents, Ai and A2, corre- 
spond to the stretching and the folding dynamics of the 
return map V c \{n) x V c i(n + 2). The spectrum that corre- 
sponds to the attractor of Fig. [3] is Ai = 0.30 ± 0.08 and 
A2 = — 0.10±0.04, with one positive Lyapunov exponent, 
and therefore chaotic behavior. 

In order to have a better picture of the way the sad- 
dle points collide with the torus we perform numerical 
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FIG. 4: We show 7 x \ f — f c \, in a log x log graph. By fitting 
the data set, we have that the average laminar time scales as 
t<x |/-/c| M , where p = 0.96. 




FIG. 5: (Color online) We show the torus (gray dots), the 
unstable focus (light gray square), the saddle points (filled 
circles), and a picture of the manifolds (diamonds) of the 
saddle points. These manifolds are depicted by iterating a 
set of initial conditions very close to one saddle point. The 
stable manifold of the saddles are indicated by dark diamonds 
(see the dark arrow) , and the unstable manifolds are indicated 
by the gray diamonds (see the gray arrow). 



analysis of the circuit equations, with parameters that 
adequately reproduce the experimental scenario. The 
circuit equations can be obtained by applying the Kir- 
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choff's laws. The resulting state equations are given 

^ : = ¥ V c2 ~ Vol) - iNR(Voi), C 2 ^f = 

^(V c i-V c2 )+i L , = -Va-Vpsmh and § = 2nf p , 

where V& and V C 2 are the voltage in the capacitors C\ 
and C2, respectively, and %l is the electric current across 
the inductor L. 

The parameters used in the numerical simulation are 
Ci = 0.1, C 2 = 1.0, L = 1/6, i = 0.575, m = -0.5, 
mi = —0.8, and _B p = 1.0. The experimental scenario 
can be reproduced if we set the perturbation amplitude 
at V = 0.234 and change the frequency in the interval 
f p = [0.1775, 0.1795]. The numerical analysis are done in 
a Poincare section at V Cl — —1.5. Plotting V C2 (n) x f p we 
obtain similar pictures as in the experimental scenario. 

In order to detect the saddle points as well as the fo- 
cus, we introduce a numerical technique, which can be 
adapted to experimental data series. The method con- 
sists in integrating the equations for several trials and 
changing the initial condition in each trial. The initial 
conditions are randomly chosen within a volume cen- 
tered at a point po m the chosen Poincare section of 
the phase space. The point po has to be close to the 
searched period- n unstable orbit, the distance was set to 
be smaller than 10^ 8 . Applying this method we were 
able to detect the unstable periodic orbits of the flow, 
which gives place to the saddle points in the Poincare sec- 



tion. The unstable orbit exists for all parameters within 
f p = [0.1775,0.1795] and V = 0.234. 

At a particular parameter f p — 0.17893, we manage 
to get a local illustration of the manifolds of the saddle 
points. For that, we evolve forward and backward a small 
set of initial condition very close to a saddle point, just a 
few returns to the section. Then, having a better picture 
of the local structure of the manifolds, we choose new 
initial conditions, which are close to the manifolds and 
iterate them. The result is shown in Fig. with the 
saddles indicated by the filled circles, coexisting with the 
pentagon-shape torus T 2 in gray dots. In this figure one 
can clearly see the distortion of the torus T 2 due the 
manifolds of the saddle points. 

Concluding, we analyze experimentally the scenario of 
a two-frequency torus breakdown by a global bifurcation. 
We show: (i) The existence of the saddle points; (ii) 
The existence of the focus point; (iii) a local picture of 
the manifolds of the saddle points; (iv) The verification 
of a power scaling law for the average laminar length, 

T0C|/-/ C |- - 96 . 
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